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Abstract 

We introduce and study an integrable boundary flow possessing 
an infinite number of conserving charges which can be thought 
of as quantum counterparts of the Ablowitz, Kaup, Newell and 
Segur Hamiltonians. We propose an exact expression for overlap 
amplitudes of the boundary state with all primary states in terms 
of solutions of certain ordinary linear differential equation. The 
boundary flow is terminated at a nontrivial infrared fixed point. 
We identify a form of whole boundary state corresponding to this 
fixed point. 
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1 Introduction 

The so called hairpin model of boundary interaction was introduced in [1] . 
This two-dimensional model of Quantum Field Theory (QFT) involves a 
two-component Bose field X(cj, r) = (X((T, r), y(cr, r)) which lives on the 
semi-infinite cylinder t = t -\- lirR, a > 0. In the bulk, o" > 0, the field X 
is a free massless field, as described by the bulk action 

M,uik = - / dr / dadX-BX, (1) 
71" Jo Jo 

with d = ^ {da — idr) and d = ^ {d^ + idr). The boundary values of this 
field, X|o-=o = i^B, Yb), are subjected to a nonlinear constraint 

(^) - (^) = , (2) 

where n is a positive parameters. A remarkable feature of the model is 
that it possesses an extended conformal symmetry with respect to certain 
W^-algebra. The generating holomorphic, Wg = Ws{cr + it), and antiholo- 
morphic, Ws = Ws{cr — it), currents of this algebra have spins s = 2, 3, 4 . . . 
and s = —2, —3, —4 . . . respectively. Among them there are spin-(ib2) cur- 
rents which are components of the stress-energy tensor: 

W2 = -dXdX -dYdY + ^d^X , (3) 
W2 = -BxBx - ByBy + ^ B'^X . 

The first nontrivial holomorphic current reads explicitly as follows 

= ^ [BYf + 2n [BxfdY + 

nlfRB'^XdY -{n^2)^dXB'^Y B^Y , (4) 

while the higher currents W4, VF5 . . . can be generated recursively from the 
singular parts of operator product expansions of the lower currents. The 
antiholomorphic currents Wg can be obtained from the corresponding holo- 
morphic one by means of the formal substitution B ^ B. 

As usual in QFT, an effect of the boundary can be described in terms 
of the boundary state which incorporates all information about boundary 
conditions [2-5]. In our case the boundary state | B )hair is a special vector in 
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the space of states TC of the two-component uncompactified scalar associated 
with the "equal time section" a = const: 

|5)haireH = ^^P®J^P , (5) 

where J-^p (-^p) is the Fock space of two-component right-moving (left- 
moving) boson with the zero-mode momentum P = (P, Q). The above 
mentioned T^-invariance of the boundary condition @ implies that the cor- 
responding boundary state obeys an infinite set of equations [3,4]: 

[Ws+lir] - Ws+i{t)]^^^ I i3)hair = . (6) 

Once preserves the conformal symmetry, the hairpin boundary condition 
can be treated as a Renormalization Group (RG) fixed point in the space of 
boundary interactions of two-component free Bose field. Broadly speaking 
any relevant boundary perturbation breaks down the scale invariance of the 
original model and introduces some RG invariant "physical scale" E^, in 
the theory. Unfortunately there is no systematic machinery for study an 
arbitrary perturbation. Therefore it makes sense to consider a particular 
class of perturbations preserving some amount of the original VF-symmetry. 
The boundary state for such models satisfies the conditions 

(I, -i,) |S)pert = (7) 

for some operator- valued functionals Ig (Is) of the original holomorphic (an- 
tiholomorphic) VF-currents of the form 

where the densities Pg+i = Ps+i[W2, W3 . . Pg+i = Ps+i[W2, W3 . . .] are 
appropriately regularized polynomials in PF-currents and their derivatives. 
Here the subscript s + 1 labels the spin of the local field Pg+i- Roughly 
speaking, the meaning of is that the boundary neither emits nor absorbs 
any amount of the combined charges Ig — Is • For this reason we shall call Is 
as a local Integrals of Motion (IM) of spin s. Notice that in the case 

P2 = W2, (9) 

where W2 is the holomorphic component of stress-energy tensor, Eq.Q man- 
ifests the invariance with respect to translations along the r-direction. 
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Let us assume that an infinite sequence of polynomials {Pg+i}, such 
that the associated IM are mutually commutative, 



[Is, Is'] = , 



(10) 



is given. It is natural to expect that a theory possessing such an infinite 
commuting set is integrable [5]. 

At the best of our knowledge the complete algebraic classification of 
infinite commuting sets of local IM for the hairpin M^-algebra has not been 
obtained yet. Nevertheless at least three nontrivial examples are known 
[6-8]. In Ref. [1] it was studied the RG boundary flow associated with the 
so called "paperclip series" of local IM. The series contains local IM with 
the odd spins s = 1, 3, 5... . The second known Abelian subalgebra [7] 
is deeply related to the Ablowitz, Kaup, Newell and Segur (AKNS) soliton 
hierarchy [9,10]. In fact the corresponding are the quantum counterparts 
of the AKNS Hamiltonians. For this reason we shall refer to this infinite 
sequence of commuting integrals as AKNS series. Among characteristic 
properties of this series is that it contains the local IM with s = 1, 2, 3 . . . 
and the first two local densities are given by Eq.® and 



One more series of the local IM containing the odd spins only is known. 
Despite admissible spins of Is for this series and for the paperclip series are 
the same, they are not equivalent. An explicit form of I3 from the third 
series can be found in [6]. 

In this article we study an integrable model associated with the AKNS 
series of local IM. The theory can be defined by adding some special bound- 
ary term of the form 



to the bulk action The potential U{X.b) turns out to be unbounded and 
pure imaginary. In spite of these somewhat pathological properties, the cor- 
responding QFT appears to be well defined and possesses many remarkable 
features. We shall refer to this theory as the Integrable Perturbed Hairpin 
(IPH) model^ 

As was pointed out in Ref. [11] a boundary state associated with inte- 
grable boundary flows with conformal bulk can be studied in a framework 

^Because of the unbound property of f/(Xs) the term "perturbed", here and bellow, 
is used in a loose sense. 



P3 = i . 



(11) 




(12) 
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of Quantum Inverse Scattering method [12-14]. To recall the basic idea 
by the example of IPH model let us assume that the corresponding local 
IM are hermitian operators and the set {1^ is "resolving", i.e., that all 
eigenspaces of the local IM are one-dimensional and mutually orthogonal 
(which seems to be the case, see Section 14. 2|) . Then, as it follows from the 
structure of space Ti © and the condition ((Tj), the IPH boundary state can 
be written as 



B 



/ d^P V5,(P) |a,P)®|«,P) , (13) 



iph 

OL 



where {|a, P)} is the orthonormalized basis of eigenvectors in each Fock 
space labeled by some index a. In a view of Eq. (|13|) . it is convenient to 
think of the boundary state in terms of the associated boundary operator. 
The natural isomorphism between T-p and JTp (the right movers are replaced 
by the left movers) makes it possible to establish one to one correspondence 
between states in !Fy' ® .^p and operators in JPp. Thus the boundary state 
I B )iph can be reinterpreted as an operator 

B= / d^P VS„(P) |Q,P)(a,P| , (14) 

which commutes with the all local IM: 

[B,I,]=0. (15) 

The structure ()13|) emphasizes an importance of the problem of simulta- 
neous diagonalization of local IM as operators acting in T^. In this con- 
nection, it is pertinent to remind that physical quantities like the boundary 
state (operator) essentially depend on some RG invariant scale E^,. For the 
IPH model it is convenient to choose this dependence in the form B = B(A), 
where the dimensionless parameter A is related with this scale by 

At the same time the associated local IM do not involve any particular en- 
ergy scale and their eigenvectors, {| a, P )}, do not depend on A. Hence the 
operators B (|14|) corresponding to different values of the "spectral" param- 
eter A commute between themselves: 

[B(A), B(A')] =0 . (17) 
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The problem of simultaneous diagonalization of commuting operator 
families is a typical problem for the Quantum Inverse Scattering method 
[12-14]. In this approach commuting families, like ]B(A) ()17() . are defined 
in terms of certain monodromy matrices associated with an auxiliary linear 
problem where A plays a role of spectral parameter. In Section|21 we briefly 
discuss the IPH model in the classical limit to clarify its relation to the 
AKNS soliton hierarchy. Later, in Section 18.21 we shall present arguments 
that the boundary operator ]B(A) in the quantum theory can be treated as 
a version of Baxter's Q-operator [15]. In particular, it satisfies the famous 
Baxter T — Q equation: 

B(A)T(A) =]B(Aq) +B(Aq"^) (q = e-^), (18) 

where the transfer-matrix T(A) is a trace of quantum 2x2 monodromy 
matrix for the auxiliary AKNS linear problem. The operator T(A) can be 
thought of as a generating function for the AKNS series of local IM. 

It should be emphasized that the paper does not contain a rigorous 
quantization procedure of the AKNS hierarchy. It is devoted to study of the 
simplest boundary amplitude 

Z = iph(5|P), (19) 

where | P ) G W is the highest vector in the Fock modulus J--p J^p cor- 
responding to an arbitrary P. A wealth of data about Z can be obtained 
through perturbative analysis in the weak coupling domain, and by looking 
into various limiting cases of the model; Sections 1511^ are devoted to these 
tasks. Using these data we propose in Section El an exact expression for the 
vacuum overlap in terms of solutions of certain ordinary differential 
equation. Only in Section |H1 examining properties of the vacuum ampli- 
tude Z, we reveal general integrable structures, like Baxter T — Q equation 
(|18p. inherent in the theory. We conclude the paper with a discussion of an 
infrared fixed point of the boundary flow. 



2 The classical IPH model 

Before going over to QFT, we will explore the classical limit of the model. 
For this purpose it is convenient to use the field 

x=(x,y) = ^. (20) 
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Indeed, if one rewrites in terms of x the parameter n appears in- front of 
the Gaussian action. This ahows one to interpret ^ as the inverse Plank 
constant. 



2.1 Integrable perturbation of the classical hairpin model 

As n — > oo, X and y become classical fields subjected by the boundary 
condition 

[e^-cos(y)]^^o = 0. (21) 

The classical equations of motion in the unperturbed hairpin model include 
the bulk equations, Ax = Ay = 0, as well as the boundary equation 

a,x-t|,=o = 0, (22) 

where t = (— tan(y), 1) is a tangent vector to the curve (|21|) . To take into 
account the zero mode, we shall consider the classical solutions x((t, r) such 
that 

x((T, r) — ^ o" as cr ^ +00 , (23) 
R 

where ^ = {(,xi (,y) is some constant vector. 

In the large n limit the quantum field W3 @ produces a classical holo- 
morphic current, W3 — > na W3, with 

W3 = 2 {dyf + 2 {dxf dy + d'^x dy - d^y dx , (24) 

and Eq.@ for s = 3 implies that the difference W3 — vanishes at the 
boundary. One can indeed check that the last condition holds in virtue of 
the classical equations of motion and the boundary constraint ()21|). 

Now let us analyze an effect of the boundary potential ()12() in the classical 
theory. We still assume the boundary constraint (|21() . so the perturbation 
modifies the classical boundary equation of motion only: 

a,x • 1 1^=0 = / , (25) 

where / = Vx^f^ • t. For an arbitrary function U(X.b) the holomorphic 
current W3 does not generate a conserving charge because 

[ws-ws]^^^ = I - -^{f{x)dax) +3f cot{y)drX + 

2{l^ + 2f)drxd^x] . (26) 

J (7=0 
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However, if we adjust the boundary potential in such a way that (|26p can 
be written in the form 



with some boundary field 62 = 02{t), then the charge 

/>oo 

g2= / da{w3iT,a)+tlJs{T,a)) -e2iT) , (28) 
Jo 

will not depend on r [5]: 

^ = 0. (29) 
dr 

This occurs for / ~ e~^^ = ^^Ji(^y\^ , i-e., for the boundary potential such that 

TlC , . , , 

U — — tan(?/) as n — > 00 , (30) 

where C is an arbitrary dimensionless constant. 

If one takes r to be the Euclidean time, Ea. H29() signifies a presence of 
nontrivial conservation charge in the perturbed theory. As a matter of fact 
the classical hairpin model with the boundary potential ()3U() possesses an 
infinite number of conserving charges Qs with s = 1, 2, 3... . They are 
generated by classical holomorphic currents Wg+i satisfying the condition 

[ws+i-Ws+i]^^Q = i^ . (31) 

For s = 1, 

W2 = -{dxf - {dyf , (32) 

and the corresponding charge qi coincides with the energy. By means of 
direct calculation it is not hard to find an explicit form of w^: 

Wi = -{dxf -f>{dyf -Q{dxf{dyf -Ad'^x{dyf + Ad'^ydxdy + 

d^xdx + d^ydy . (33) 

It is useful to keep in mind simple ambiguities in a choice of Ws+i satisfying 
Ea. (|ST|) . First, these currents are defined up to total derivatives, Ws+i — > 
Ws+i + dg with dg = 0. Second, the equation (|ST|) does not fix an overall 
multiplicative normalization of Ws+i- 
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2.2 Local IM in the AKNS soliton hierarchy 

Now we describe an effective way to generate all the densities Wg+i (|31jl up 
to the above mentioned ambiguities. At this step we need to introduce the 
fields 

il) = \ {dy + \dx) exp{2iyR) , ip* = i {dy -idx) exp(-2iyi?) , (34) 

where yn in the exponential stands for the holomorphic part of the harmonic 
field y: 

y{T,cr) = yR{(T + \t) +yL{a -it) . (35) 

Hence il) and ip* are non locally expressed in terms of the fundamental field 
X. They are clearly holomorphic {dtp = dip* = 0) and quasiperiodic fields: 

iP{t + 27tR) = e-^'^'^y V(t) , '^*{t + 2ttR) = e^^'«^ tp^r) , (36) 

where S^y is the second component of the vector ^ in Ea. (|23() . 

The fields (|34j) are remarkable in many extents. First of all, the densities 
Ws in Eas. (PH) . are nicely expressed in terms of ip and tp*: 

W2 = 'P'P* , 

W3 = i {rd^ - ipOr ) , (37) 
W4 = -\ {tP*d^ip + ipd^ip*) -{ipip*f . 

More importantly, ip and tp* generate a closed Poisson subalgebra in the 
space of classical fields. To introduce the Hamiltonian picture here, we 
will interpret the world-sheet coordinate a as the Euclidean time. Then 
the classical bulk action for the fundamental field x defines a canonical 
Hamiltonian structure which implies the following set of Poisson brackets 
for the holomorphic components: 

{xr{t) , xr{t')} = {yR{T) , yR{T')} = -\ e{T - t') , 

{xr{t) , yR{T')} = 0, (38) 

with 

e(r) = 21 + 1 for 2t:R I <t < 2ttR (/ + !); leZ . 
Using Eqs.(|2Hl) it is easy to show that 

{V'(t), V'(r')} = e(T-T') V(t)V'(t'), 

{ V*(r) , ) } = e(r - r') ^(t)^ (r') , (39) 

{ r{r) , V'(t') }=S'{t- t') - e(r - r') r{r)iP{T') . 
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The Poisson bracket algebra H39|) is well known to describe the second Hamil- 
tonian structure for the AKNS soliton hierarchy [16, 17] and Eas. ()34() can 
be interpreted as a transform to the Dorboux variables ()38() for this Hamil- 
tonian structure. 

The fields (|37|) are local densities of the first AKNS Hamiltonians: 

^(ciass) ^ ji-. dTu;,+i(r) . (40) 



Jo 

In particular j^'^^^^^ coincides with the Hamiltonian of the famous non-linear 
Schrodinger equation [9] ^. There is an infinite number of IM in the form 
H40() which form a commutative Poisson bracket subalgebra 

|^(class)^ j(class)^^Q^ (41) 

It is remarkable that all AKNS local densities Ws (|4(H) satisfy Eq. (|31j) . This 
statement is essentially known in the literature in the context of Backlund 
transformation (see, e.g.. Appendix B in [18]), even though it is not formu- 
lated in the language of the perturbed boundary theory. 
It is well known (see, e.g., [17]) that the AKNS flows, 

g = { 4^1-) , ^'(ti, t2 ...)}, ^ = { 4^1^-) , r{tut2 ...)}, (42) 

describe isospectral deformations of the first order differential operator 

C = -dr-'4 H + Tp* E + iP F , (43) 

where tp, ip* are the quasiperiodic fields on the segment < r < lirR 
and E, F and H are the generators of the Lie algebra sl{2), 

[H,E] = 2E, [H,F] = -2F, [E,F]=H. (44) 

More precise, if we define the (2j + 1) x {2j + 1) monodromy matrices 
Mj(A), j = 0, 1/2, 1, 3/2..., corresponding to the (2j + l)-dimensional 
representation ttj of (|44j) . as 

(Xi(r + 27ri?) , . . . X2j+i{r + 2vri?) ) = (xi(r) , . . . X2j+i(.r) ) M,(A) , (45) 



^The classical IM in Eg. 1141)1 are normalized in accordance with the convention from 
the book [17]. If {xr, i/r) in H84|l are real functions of the real variable r, then tp and tp* 
is a complex conjugated pair and the corresponding non-linear Schrodinger equation 1421 1 
is in the repulsive regime. 
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where Xki^) (2j + 1) linear independent solutions to the auxiliary linear 
problem 

7rj[C]x = 0, (46) 

then the transfer-matrixes 

r,-(A) = Tr^^. [e~^^'^y^ M,(A) ] (47) 

are involutive (with respect the Poison structure (jSHJ) of the AKNS 
flows: 

{T,(A), 4^'^^^)} = (48) 

and 

{T,(A),T,v(^)} = 0. (49) 

As the matter of fact the transfer-matrices Tj{X) are not independent for 
different j and can be algebraically expressed in terms of T(A) = T'i(A) 

2 

corresponding to the fundamental representation of sl{2). The latter can be 
thought of as a generating function of the local IM (|4fl|) as it expands in the 
A — > oo asymptotic series [17]: 

r(A) = 2 cosh (27ri/(A)) , (50) 

with 

oo J-. 

i -A - + 2^ ^ jcciass) ^ ^ ^ ^ _ (51^ 

s=l 



3 Semiclassical quantization 

Here we study the semiclassical behavior of the boundary amplitude (|19() 
using the path integral approach. For the sake of discussion it is convenient 
to consider the conformal map of the semi- infinite cylinder, r = t+IttR, a > 
0, to the disk of radius R: 

- = e-(-+'-)/^ , ^ = e-('^-i^)/^ (52) 
R R 

Then the overlap pert ( -B | P ) with the Fock vacuum | P ) relates to the un- 
normalized one-point function of associated primary field inserted at the 
center of the disk, 

{^^■^{0,0)),,^, = Ry'-^'/\^,,{B\P), (53) 
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and 

pei.t(i? I ) is the disk partition function"^. The one-point function 
H53() can be represented in terms of the functional integral as follows, 

(eiP-^(0,0) )^.^^ = J VXVY ei™«^(0,0) e-^I^l-^p-'I^B] , (54) 

where P = (P, Q), and the integration variables X{z, z), Y{z, z) are assumed 
to obey the constraint © at the boundary \z\ = R. Notice that at the 
classical limit the effect of exponential insertion in H54|) can be accounted 
for by imposing the asymptotic condition H23() in the cylindrical frame. For 
this reason the vector ^ in H23|) is proportional to P: 

^ = ^ • (55) 

When P and Q are pure imaginary, some insight can be gained by making 
a shift of integration variables, 

X^X + iPlog^, (56) 
in the functional integral (|54j) . which brings it to the form 

(eiP-^(0,0)>^j^j^ = ii-P'/2 j VXVY e-^buik[x]-^bow[Xfl] ^ ^57) 
where the boundary action is given by 

-/ ^ {PXB + QYB + iRU{XB,YB)){z) . (58) 



'bound 



In this section we evaluate leading semiclassical contribution to the overlap 
H19() by direct calculation of the functional integral H57() in the saddle-point 
approximation. This will give some intuition about its structure. 

3.1 "Light" vertex insertion 

Making an attempt to calculate ()54() we immediately encounter a problem; 
The potential (|3n|) is unbounded as y — > it^. In fact, in order to ensure a 
convergence of the functional integral we must place the overall pure imag- 
inary constant C in (|8n|) ^: 

C = iA . (59) 

^We always assume that the exponential field in (1531 is defined according to the usual 
CFT conventions [19]. The factor B}^^'^ appears due to the conformal anomaly. 

*It deserves to draw reader's attention that in I59II we use the same notation A as for 
the spectral parameter in Eg. 1431 . The reason is discussed at the end of this section. 
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For real positive A, the path integral ()57() seems to be well defined, but the 
pure imaginary boundary potential U(X.b) will certainly break down the 
unitarity of the underline QFT in the Minkowski coordinates {a, it). 

Let us assume that the parameters P, Q and A are sufficiently small, 
i.e., the boundary action (|58|) has no appreciable efi'ect on the saddle-point 
configuration. To be precise we write 

(P, Q) = ^ (P, q) (60) 

and assume that 

p, q and nA ~ 1 as n ^ oo . (61) 

Notice that these conditions are somewhat different then those considered in 
the previous section, where the limit n — > oo was taken under the assumption 
that = p/n, = q/n and A are fixed. 

With the condition (|61j) the action is minimized by the trivial classical 
solutions X(z, z) = Xq, where Xq = (Xo,lo) is an arbitrary point on the 
hairpin Q. Therefore 

Zciass = / dA^(Xo) e2i(P^o+<?yo)/v^ e-RUi^o) ^ (62) 

./hairpin 

where the integration measure dA4(Xo) is determined by integrating out 
fluctuations around the classical solution in the Gaussian approximation. 
Of course there is no need of actually evaluating this Gaussian functional 
integral to figure out the answer. If one writes 'X.{z,z) = Xq + 5X(z, z), and 
splits the fluctuational part (5X into the components normal and tangent to 
the hairpin at the point Xq, the components are to satisfy the Dirichlet and 
the Neumann boundary conditions, respectively. Therefore one just has to 
take the product of the known (see, e.g., Appendix to [20]) disk partition 
functions with these two boundary conditions. As the result, dA^(Xo) = 

ll d^(Xo), where d^(Xo) = v^(dXo)2 + {dYoy is the length measure of the 
hairpin, and g_D = 2~4 is the "boundary degeneracy" [21] associated with 
the Dirichlet conformal boundary condition^. 

Certainly the above analysis ignores the presence of (unbounded) bound- 
ary potential. Its proper account leads to the one loop renormalization of 

^The definition is as follows; go ~ {Bd\P), where \Bd) is the boundary state of 
uncompactified boson X with the Dirichlet boundary condition Xb = 0, and the primary 
states I P> are delta- normalized, (P | P' ) = 5{P - P'). 
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the couphng constant A. Indeed, the form of the boundary potential in the 
classical limit suggests that 

U{Xb) ^ ±i — e ^ as ^ -oo , (63) 
it 

where the sign factor is dictated by the choice of hairpin asymptote Yb = 
=p^^==. The anomalous dimension of the boundary operator e ^ with 
respect to the hairpin stress-energy tensor ^ is given by 

dpert = . (64) 

If one imposes now the normalization condition 

^ {XB)disk] =0 (65) 



E,=0 



on the expectation value of renormalized boundary field Xb, then the renor- 
malized coupling A is related to the RG invariant scale Ej, as in Eq. (|16|) . 

Although Ea. H16|) is a result of the semiclassical consideration, in the next 
section we shall bring forward arguments which show that it is perturbatively 
exact, i.e., a renormalization scheme exists in which it is exact to all orders 
in the ^-expansion. 

All the above add up to the following integral representation for the 
semiclassical boundary amplitude: 



^class 



gl ^ (f r^A)-'^ j\ g (cos(y))^'^-^2'^^ e'"^*^'^(^) , (66) 



where the additional factor {'^r'^X)^^^ containing an arbitrary constant r 
is introduced to ensure the normalization condition H65|) . This factor can 
be also understood as an artifact of an additive counterterm for the bare 
boundary field Xb- The integral H66|) is expressed (see Eq.l3, Section 6.11.2 
in Ref. [23]) in terms of Kummer's solution U{a,b;z) of the confluent hy- 
pergeometric equation [24]: 



^ciass = gl r-''^ ^^""^^ e-"^ [/( 1 _ + 1 _ 2ip- 2nX ) . (67) 
T{^-q + ip) 

Here we assume that A > and choose the principal brunch of the multi- 
valued function U (a, b; z) such that for positive z 

U{a, b; z) z'" as z^ +oo . (68) 
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It is also instructive to rewrite the semiclassical boundary amplitude (|67|) in 
terms of the hypergeometric series M{a, b, z) = 1 + ^ + ^^ij'j^l^ z'^ + . . . : 

^ciass = 4;L (p^ 9) Pc^2s (p^ 9 1 A) + 4tL (p^ <?) ^iiti (p> i a) , (69) 

where 

Z^-Uva) - ^'r-^^ ^ ^^^'^^ (70) 

Z^,^,{p,q) - gnr + + (70) 

zitliP^q) = glr-^'P ^i2nXyP cosh {nip + iq))Ti-2ip), 

and 

^iili (P, 9 I A) = M ( i ^ ip + g, 1 ^ Tip; 2nX) . (71) 

In this form the nature of singular behavior as A —> is more explicit. Let 
us make a (trivial) observation that the poles of ■^ciass(^'' l) variable p 

in the first term in (|69() are canceled by the poles in the higher terms of the 
confluent hypergeometric series Pj:^^ in the second term, and vice versa, it 
makes the full partition function an entire function of parameters p and q. 



3.2 "Heavy" vertex insertion 

Although the above result 1)69^ was derived under the assumption that P,Q 
and A are small, it needs little fixing to become valid for much larger values 
of these parameters. When {P,Q) = (CziC?/) become as large as ^/n, 

^x, aiid A ~ 1 as n — > oo , (72) 

the vertex insertion in (|54|) and the boundary potential is "heavy", i.e., it 
must be treated as a part of the action, and they affect the saddle-point 
analysis. The saddle-point configuration(s) is still a constant field, but now 
Xq is not an arbitrary point on the curve but has to extremize the 
boundary action (|58|) 

[Xo] = -iP-Xo + i2C/(Xo) . (73) 

There are two solutions of the saddle-point equation which are real valued 
for real values of the parameter 

uo = ^il-2\iy-X^ . (74) 
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One of them corresponds to the minimum, and another to the maximum of 
i =g4ound[Xo]- The saddle-point action produces nothing else but the p, g — > 
oo asymptotic form of the expression H69() - after all, this asymptotic of the 
integral is controlled by the same saddle points. One can observe that 
if one splits the constant-mode integration into two parts, as was suggested 
in H69|) . the parts receive contributions from different saddle points - one 
from the "minimum" and one from the "maximum" one. What makes the 
difference at {P, Q) ~ ^/n is the proper treatment of non-constant modes. 
One writes 

X(z, z) = X* + t* 6Xt{z, z) + n=K 6Xn{z, z) , (75) 

where X,,, is the position of the saddle point on the hairpin, and t,,, and n,,, 
are unit vectors tangent and normal to the hairpin at this point. Then for 
small 6Xt the boundary constraint Q reads 

-^0 SX? o 
6X^ = -e^ -^+0{6Xf), (76) 

and, up to higher-order terms, the boundary action l|58j) can be written as 

Ab = Ab[X, ] t z.0 ^ ^ SX? , (77) 

with the coefficient z/q > is given by (|74l) ; then the sign minus (plus) in (|77|) 
applies to the "minimum" ("maximum") saddle point. Thus, while to the 
leading approximation the normal component 6Xn still can be treated with 
the Dirichlet boundary condition, the "boundary mass" term in (|77j) has to 
be taken into account in evaluating the contribution from 6Xt. Using the 
well known boundary amplitude of the free field with quadratic boundary 
interaction [22] , one finds that Eq. H69|) would apply to the case of (P, Q) ~ 
-y/n as well if one puts corresponding additional factors in the two terms in 
(|69|l . i.e., replaces Z^'^^^{p,q) there by 



41(^> Q) = 4lL(f P , f Q) ri ± i J^-^X-4X- ; (78) 



of course in this case one can use the p, g — > oo asymptotic forms of the 
factors (f7fl|) and (f7T|) . More explicitly, 

^c,a. - 4= ( ^-^^ e-- + ^-^^^ e-^ \ , (79) 
V2^ V V2i^ V^2i^ J 
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with 



A 



dS± 



(80) 



dX 



Here the two terms in H79() correspond to the different saddle points of the 
boundary potential (|73|) . The structure (|79|) is common for the semiclassical 
form of Baxter's Q-function (see Ref. [25]). It is instructive to note in 
this connection that the classical transfer-matrix (A) corresponding to 

simplest ("vacuum") choice x/j = ^ (cr+ir) in Eas. (|51)) . (jl^ . can be written 
in the form ()5U|) : T(™'^)(A) = 2 cosh(27rt'o), where uq is the same function 
(|74j) as in Eq. (|8()|) . One can easily recognize in these equations the classical 
counterpart of the Baxter relation (|18() . 

4 Integrability of the IPH model 

In this Section we would like to argue that the quantum IPH model is inte- 
grable. Our first step toward the consistent quantum theory is the quanti- 
zation of the local IM of AKNS soliton hierarchy. 

4.1 Quantum local IM 

Let us first recall structure of the hairpin VF-algebra [1] ^. Its generating 
currents Ws{u) (n = cr + ir) have spins s = 2, 3, 4 . . . , and can be charac- 
terized by the condition that they commute with two "screening operators 
(charges)", i.e., 



and Xf(, Yfi in the exponential stand for the holomorphic part of the cor- 
responding local fields. The integration in Ea. (|81|) is over a small contour 
around the point n; vanishing of the integral (|81|) implies that the singu- 
lar part of operator product expansion of Ws{u)V±{v) is a total derivative 
dy{. . .). This condition fixes Ws{u) uniquely up to normalization and adding 
derivatives and composites of the lower-spin Il/'-currents. For instance, the 

®This PF-algebra was known for a long while (see [6,26,27]). 




(81) 



where 




(82) 
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first holomorphic current beyond ^ can be written as @ where the ambi- 
guity in adding a term proportional dW2 is fixed by demanding that @ is a 
conformal primary. The higher currents W4, W5, . . . can be found either by 
a direct computation of the operator product expansions with the screening 
exponentials (|82|). or recursively, by studying the singular parts of the op- 
erator product expansions of the lower currents, starting with W3{u)W3{v) 
and continuing upward. Thus, the product W3{u)W3{v) contains singular 
term ~ {u — v)~'^ which involves, besides the derivatives d'^W2 and the com- 
posite operator VFf , the new current W4. Further operator products with 
W4 define higher VF's, etc. In this sense the VK-algebra is generated by the 
two basic currents W2 and W3. 

As the matter of fact, there is the third, the most effective way to gener- 
ate all VF-currents. It based on the observation that the screening operators 
associated with the exponentials 1)82(1 commute with the parafermionic cur- 
rents 



^{u) =i(v'^rr2 ay-hi^/^5X) ev^^^^"^ (83) 
^*{u) =i(^/^rr2 ay-iVr^5X)e"^^^^''^ , 
in the same sense as the W^-currents do, i.e., 

I dv^{u)V±{v) = Q, (f dv ^*{u)V±{v) = 0. (84) 

J u Ju 

The fields (|83|) extend the notion of the parafermions of [28] to non- 
integer k = —n — 2. Clearly they can be also treated as a quantum version 
of the classical nonlocal currents (|34|) . In spite the fact that the fields (|84|) 
are not local, both ^ and ^* are local with respect to the exponentials 
(|82|) . hence the integration contour in (|84|) - a small contour around n ~ is 
indeed a closed one. It follows from (|84() that all the fields generated by the 
operator product expansion of ^'(n)^'*(ti) satisfy Eqs. ((HT|) [6,27]. Thus we 
have 

^{u)^*{v) = (n - v)-^^ { + § {W2{u) + W2{v)) - 

{Ws{u) + Ws{v)) +...}, (85) 

where W2 and W3 are the same as in (j^J and @ , and the higher-order terms 
involve the higher-spin Pi^-currents. 

After the brief review of the hairpin extended conformal symmetry, we 
turn to the description of the quantum AKNS local integrals which consti- 
tute the Abelian subalgebra of the VF-algebra. The main idea is based on 
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examination of the asymptotic form of classical boundary potential 1)63(1 : 
It suggests to define the quantum local IM Ig as elements of the hairpin 
VF-algebra in the form ^ commuting with the additional screening charge, 
i.e., 

(f dv Ps+iiu)Vo{v) = du{ ...) , (86) 

where 

Vo = e"^ . (87) 

Remarkably the condition (|86|) indeed defines, up to an overall multi- 
plicative normalization, a local integral of motion Ig for each s = 1, 2, 3 . . .. 
In particular, one can show that P2 = W2, P3 = iVFs and 

Pi = n (dX)'^ + 6n {dXf{dYf + (5n + 4) (dY)^ + (88) 
6(n + l)^/^^'^X{^Yf{n^ + 3n + 1) {d'^Xf + (n^ + 4n + 2) {d'^Yf, 

where we disregard all terms which are total derivatives and do not con- 
tribute to the I3. Also it is straightforward to check that the such defined 
Ii, I2 and I3 are mutually commute. One may also note that in the classical 
limit P4 (|55|) turns to be the classical local density (— ^^4) (1^5]) up to a total 
derivative: P4 ^ ( — W4 + d{. . .)) as n — > 00. 

It seems likely that the formal proof of existence and uniqueness solution 
of ((86|) for s = 5, 6... can be obtained along the line of Refs. [29-31]. 
Currently it is not available. Nevertheless we take the foregoing facts as a 
strong indication that an infinite sequence of commuting IM { in the 
form Q exists, whose first local densities are given by Eqs. ©, (|TT|) and 
(|88() . It is expected that the operators Ig become the AKNS Hamiltonians 
(|4()|) in the classical limit. 

4.2 Diagonalization of local IM 

It seems natural to take an attitude that the boundary state of quantum 
IPH model satisfies the integrability condition ^ for all the local IM from 
quantum AKNS series. Thus we run into the problem of simultaneous diag- 
onalization of the AKNS local IM in the Fock space J-p . 

Let us consider the normal mode expansion of the quantum holomorphic 
current 9X: 

dX = j^ (ip + ^Xfce^^ ). (89) 
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Here P = {P,Q) is the zero-mode momentum and = {Xk,Yk) are the 
oscillatory modes. The canonical quantization procedure of the Poisson 
bracket algebra H38|) leads to the following set of commutation relations, 

[X^Xr] = |5^'^ W- (90) 

In Appendix A we present explicit expressions for the first local IM in terms 
of the oscillatory modes Xj^. As it follows from Eas. H211|) . (|212() . are well 
defined operators acting in the Fock space J^p. The later is the highest 
weight module over the Heisenberg algebra (|9U() with the highest vector 
I 0, P) ("vacuum") defined by the equation 

|0, P) =0, A; = 1,2,3... . (91) 

The zero- mode momenta, P and Q, acts in the given Fock space as c- 
numbers. The space J-j> naturally splits into the sum of finite dimensional 
"level subspaces" 

= eT=o^P ; = ^ ^P ■ (92) 

Since the grading operator L essentially coincides with Ii, 

Ii = i?-i (L+^-3^) , (93) 

{(■) 

all the local IM act invariantly in the level subspaces J-^ . Therefore diago- 
nalization of Ig in a given level subspace reduces to a finite algebraic problem 
which however rapidly becomes very complex for higher levels. Notice that 
for n > and real P, Q, the operators are hermitian with respect to 
the canonical conjugation {X^)'^ = X'^f^, so their spectra are real and their 
eigenvectors are orthogonal to each other. Here we list only the vacuum 
eigenvalues for the first few I^: 

I, I 0, P ) = /i™^) (P, g) I 0, P ) , (94) 



where 



.(vac) ^ El>9l_l 

1 4 ^ 4 12 ' 

(5n + 4)(g^-g^ + X)^(!i±^^. 
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An explicit form of the eigenvalues and the corresponding eigenvectors in 
J^p is somewhat cumbersome even for small i. For this reason we do not 
present it here. However, it suggests that even the first three local IM resolve 
all degeneracies in J^p, i.e., they orthonormalized eigenvectors {| a, P)}, 
which are labeled by some index a, form a basis in the Fock space. This 
is in turn leads to the structure Notice that the amplitude Bo(P) in 

Ea.()13|) corresponding to the normalized highest vector | 0, P) G J^p, 

(P', 0|0, P) =(5(P'-P) , (96) 

coincides with the complex conjugated overlap Z* ()19() . 



5 Dual form of the IPH model 

Here we discuss an alternative description of the IPH model and explore 
the short distance behavior of the boundary amplitude H19() beyond the 
semiclassical approximation. 



5.1 Dual Hamiltonian 

In the recent work [32] it was suggested that the hairpin model admits equiv- 
alent "dual" description. The dual model involves a two-component Bose 
field {X{a, r), y (a, r)) (where Y is interpreted as the T-dual of Y on the 
semi-infinite cylinder, which has the free-field dynamics in the bulk, and 
obeys no constraint at the boundary a = 0; instead it interacts with an 
additional "boundary" degree of freedom. In the dual description it is con- 
venient (but not necessary) to use the Hamiltonian picture where the cyclic 
coordinate r = r + 2ttR is treated as the Matsubara (compact Euclidean) 
time T. In this picture the boundary amplitude Z^~\P,Q) = hair(-B|P) 
admits the dual representation as the trace 

Z^~\P, Q) = Tr^ [ Q-^^RH^^ir 1 ^ (-97) 

taken over the space Ti = 'H-^y ® where 'H-^y is the space of states 
of the two-component boson (^X{a),Y {a)) on the half-line u > (with no 
constraint at C7 = 0) and C'^ is the two-dimensional space representing the 

^The T-dual of free massless field is defined as usual, through the relations 
Y = Yr{g + ir) + Yl{(J--it), Y = Yr[(j + ir) - Yl[(J - it) . 
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new boundary degree of freedom. The dual Hamiltonian in 1)9 7|) consists of 
the bulk and the boundary parts, 

-^hair = -H^bulk - 2^7? - Qo'3 + (98) 

The bulk part, i^buikj is a "free" Hamiltonian corresponding to the bulk 
action The boundary term describes coupling of the boundary values 
{XB-,yB) = {X,Y)\(j=Q of the field operators to the additional boundary 
degree of freedom represented by {a± and 173 are the Pauli matrices 
acting in C^). We also include in Ea. (|98|l an extra parameter /i_ > 0. In 
the unperturbed hairpin model /i_ can be always eliminated by shifting of 
the field X, but this coupling will be of considerable importance for the IPH 
model. 

We emphasize the similarity of the boundary vertex operators in Eq. (|98)) 
to the screening charges ()82p determining the extended conformal symmetry 
of the hairpin model. To construct the dual Hamiltonian for the perturbed 
theory, it is crucial to observe that the local IM Ig defined through the 
condition ()86|). also commute with the additional "dual" screening charge: 

idvPs+i{u)V^^'''''\v) = du{...), (99) 

J u 

where 

^(dual) ^ ^-2v^X« _ 

For s = 1, 2, 3 Ea. H99|) was checked by the direct calculation using the 
explicit forms of . We take this as a strong indication that it holds for all 
the local IM from AKNS series. This in turn suggests a possibility of dual 
description of the IPH model by means of the Hamiltonian i^iph : 7i ^ TC, 

i^iph = i^hair + 5] e-V^^s , (101) 

where S is some 2x2 matrix acting in the C^-component of the Hilbert 
space TL. In what follows we will argue that the matrix S has the diagonal 
form, 

S = /i+e-'^'^s ^ (102) 

where /i+ is an arbitrary constant. Notice that the Hamiltonian -ffiph ()101|) . 
H102|) is not hermitian. This is consistent with our previous observation from 
Sectional that the IPH model, once considered in the Minkowski space, is a 
nonunitary QFT. 
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5.2 Short distance expansion 

To justify the conjectured form of dual Hamiltonian we shall study the short 
distance behavior of the partition function corresponding to H1U1|1 : 



Tr^ 



(103) 



Qualitatively, one may expect that if SJm P is not too small, the limit i? — > 
of the partition function Z is controlled by either the hairpin boundary oper- 

y/n -tr _|_- -v/n + S -y ____ 

ators e~2~ ^ 2 b (jggj) or the Liouville boundary vertex H101|) depending 
on the sign 'QmP, with some crossover at small 9m P. 

More precisely, let us assume first that QmP < 0. Then we can treat the 
second term in (jlOlj) as a perturbation and in the leading approximation 

Z\r^o^ Z^'\P,Q) (9mP<0), (104) 

where Z^~\P,Q) = hair(-B|P) is the hairpin boundary amplitude (|97|1 
found in the paper [1]: 

nJ r(i+g + ip)r(i-g + ip) ^ ^ 

Here and bellow we often use the notations 

p=lVnP, q=\y/^^T2Q . (106) 

Notice that in the semiclassical case n ^ 1 the parameters p, q (|l()6j) are 
the same as in ()60|) ; One should not distinguish between n and n + 2 at the 
perturbative order discussed in SectionlHI 

Now we turns to the case QmP > 0, where the leading short distance 
behavior is controlled by the Liouville vertex e~^^B _ jj^ fact, because of the 
matrix factor S (\W2\i , there are two noninteracting boundary Liouville field 
theories with the boundary "cosmological constants" fi+e~~ and ;U+e~2~. 
Using the result from the works [33,34] one has, 

Z\r^o^ Z^+\P,Q) (9mP>0), (107) 

with 

ZM = g|, 1^ r( - 2ip )r(l - 2i £ ) (^^[^^) ( e'^'^^'^^ + e-^--^ ) . 

Here the last factor comes from the trace over C^-component of the Hilbert 
space Ti. 
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Since Z^^\P, Q) vanishes in the limit i? ^ if P is taken in the "wrong" 
half-plane (note the factors R~'- and W'^^r- in (|Tn^ . (|Tn71) ). this in 
turn suggests that the overall -R — > asymptotic of the partition function is 
correctly expressed by the sum 

Z\ii^o^ Z^-\P,Q) + Z'^+\P,Q). (108) 

Remarkably that Eq. (|108j ) is in a perfect agreement with our semiclassical 
result (jini), ((ZHl) provided 



//+ r(i - n) 

and 



r 



RY , (109) 



= i2E.r^' , (110) 

Ibvr^n 

where -E,,, is the RG invariant "physical scale" in the IPH model ()16() . It 
gives a strong support to the conjectured Hamiltonian H101|) . 

What can be said about corrections to the leading asymptotic? Again 
we assume that P < and consider the perturbative effect of the second 
term in H101|) to the hairpin boundary amplitude 1)104(1 . In the unperturbed 
hairpin theory the parameter /i_ is a dimensionless constant. Let us elimi- 
nate fi- from the hairpin Hamiltonian by shifting of the field X. Then the 
coupling ^_|_ in front of the Liouville term in Eas. ()101() . ()102() is replaced by 
Since the anomalous dimension of the boundary operator q~^^b 
with respect to the hairpin stress-energy tensor Q is given by 

dduai = -n - 1 , (111) 

then ~ E^'^'^, which is in an agreement with Ea. (|110|) . Hence we 

deduce that the perturbative corrections to the leading asymptotic (|104|) 
should be in a form of power series expansion of the dimensionless parameter 
with 

K = E^R . (112) 

The case QmP > can be analyzed similarly and one comes to the same 
conclusion about the form of perturbative corrections to the leading asymp- 
totic (tWI) . 

The power series expansion in k""*"^ becomes invisible in the limit n — > 
oo, while the consideration from Section|Sl suggests that the short distance 
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expansion should be in a form of power series of the dimensionless parameter 



n+2 



A ~ K " ()16() . This structure is neatly captured by the form 

Z = Z^-\P,Q) F^-\p,Q\k) + Z^+\P,Q) F^+\P,Q\k), (113) 
where the functions F^^^ {P,Q\k) are double power series in integer powers 

_i_o n+2 

of k""*" and n " : 

oo 

fW(p,q|k) = J2 fS^iP^Q) , (114) 

i,j=0 

with /o,o = 1- The exact splitting into two terms in (|113j) is not easy to 
justify on general grounds, but can be supported by the following arguments. 

First, recall that the such splitting in the semiclassical expression (|69)) 
corresponds to isolating contributions from two saddle points. More im- 
portantly, the full expression has to take care of the poles of the factors 
Z^^\P, Q) - one should expect that the partition function ()113() is an entire 
function of P. The expression ()113() is the simplest form fit for this job. 
The poles of Z^~^^ {P, Q) are located at the points in the lower half of the 
P-plane where i^/nP or i take non-negative integer values. At these 

■p(n + 2) 

points the factor k in Z^^' {P,Q) "resonates" with certain terms of 

the expansion (|114|) in the first term in (|113|) . Therefore, the poles of the 
factor Z^^^ (P, Q) in the second term can (and must) be cancelled by poles in 
appropriate terms of the expansion of the first term. For the "perturbative" 
poles at i^/nP = 0, 1, 2. . . this mechanism is evident in the semiclassical 
expression ()69() . The form similar to ()113|) , together with this mechanism of 
the pole cancellation, was previously observed in the boundary sinh-Gordon 
model [35] and in the paperclip theory [1]. 

With the dual representation (|lfllj) . (11031) it is possible to make some 
quantitative predictions about the coefficients in expansions H114|) . In par- 
ticular, using the approach outlined in Ref. [36] (see also [32, 37]), it is 
possible to obtain Coulomb gas integral representations for f^Q{P,Q) at 
the poles P = i m = 0, 1, 2 ... of the first term in the sum H113|) . Un- 
fortunately, for general m and i the such integral formulas appear to be 
too cumbersome for any practical use. The sole exception is the case with 
m = 0, i = 1: 

(_), _ 2»r(l-n) J{q) 
Ji,o |p=o 4^2^ 2cos(7rg) ' ^ ' 
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where 



ei-(9-f ) + e-'^^^-t) Ji(-g) + ei-(9+f ) J2(g) + (116) 



and 



Ji 
J2 



2sin(V; 



0<w<v<u<'2tz 



< u < KJ < Ti < 27r 



< u < < Ti; < 27r 



-n-l 



4sin(2i^) sin(^) 



2sin(^i^) 4sm(2ifi^) sin(Hz:^) 



— n— 1 



2 sin (^) 4 sin {'^) sin 



^ig(?)-n) 



^iij(i)-n) 



^iq(v-u) 



These integrals should be understood in a sense of analytical continuation 
in n from the domain of convergence (— 1 < Sften < 0). The phases in (|116p 
are also worthy of notice. They are generated by the additional boundary 
degree of freedom in (|98|) and make the sum of integrals H116|) well defined. 
Were these phases absent J{q) would depend on an auxiliary initial point 
of integration for the variables u, v, w. Here this point is chosen to be zero. 
For the given phase factors the sum of integrals 1)116^ can be expressed in 
terms of the generalized hypergeometric function at unity to yield 



(_) 2"r2(-n)r(i + <7) 

/l,0 \p=o ~ 



3i^2(|+g, -n, -n- 1, \ +q-n; l).(117) 

n 



6 The IPH model for n ^ 



The dual description of the IPH model is especially useful in the strong 
coupling [n 0) regime. In particular for n = the theory simplifies 
drastically and can be explored with full details. 



6.1 Boundary amplitude 

Here we consider the boundary amplitude (|19() in the n ^ limit, assuming 
the parameters p, q (|lfl6|) are fixed. For n = the field X formally decouples 
in the dual Hamiltonian H98|) . H1U1|) . To be more precise the n — > limit is, 
in fact, the classical limit for X. After the field redefinition, 

(j) = ^X (118) 



26 



this becomes particularly striking. Since quantum fluctuations of </> are sup- 
pressed as n — > we may apply the saddle point approximation to account 
their contribution. The theory is trivial in the bulk and, in consequence, the 
saddle point achieves at some constant classical field configuration which we 
shall denote as (pQ. Thus, in the limit n — > we find 

FiphU^o = fl'buik - ^ + M+e-'^° + (119) 



/ids + /i_e 2 



: Vn + 2 v> ■ Vn + 2 -y 

(T+ e' + (T_ e ' 



with 



/i = -1 + 7rn;U+e-'^«) . (120) 



One can observe that H119|) is a Hamiltonian of the one-channel anisotropic 
Kondo model [38] and n = corresponds to the so called Toulouse point 
[39]*^. It is well known that the ground state energy of the Kondo model 
suffers from a specific "free-fermion" divergence at the Toulouse point. For 
this reason we prefer to keep n as an ultraviolet regulator in Ea. (|119j) rather 
than set it to zero. Notice also that h (|12U|) plays the role of an external 
magnetic field coupled with the impurity spin. The partition function of 
Kondo model is well known (see, e.g., Ref. [40]), 



^ Kondo 



2vr gD e 



r(i + 2iRh + TTRfj.1 e'f'o) r(i - 2iRh + 7rRfj.l_ e't'o) 



(121) 



where 



UesiM = -^(l}o + ni^+ e-^o + f^e^°+nEo. (122) 
TTri 

The first two terms in the effective potential Ues came from the constant 
((/>o-depended) terms in the Hamiltonian (|119|) . The third term is the above 
mentioned free-fermion divergence and Eq is a non-universal constant with 
the dimension of energy. The saddle point configuration (pQ corresponds to 
the minima of the effective potential Ugg, which picks at 

ip 

sinh(0*) = — {(po = — nlogr -|- (/)=„) , (123) 

K 

where we use the notations r and k from Eqs. (|l()9j) . (|ll()j) . (|112)) with n ^ 0. 
Now we should take into account an effect of Gaussian fluctuations around 



*In the conventional Kondo model the field Y emerges from the bosonization of electron 
degrees of freedom and is assumed to be compactified. 
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the classical saddle-point configuration. We can write (/)(cj, t) = (/>o+50(o", r). 
Then expanding the effective potential near the minima, we obtain the 
boundary potential in the Gaussian approximation, 

1 



UMb) = const + ^ y/i^^ {S<Pb? + 0{i6(l)Bf) . (124) 

Finally we can use the result from Ref. [22] to find the contribution of the 
Gaussian field 64> into the partition function. Combining all these ingredients 
together one arrives to the following result 



ZL^o = gn r 



2 ^-2ip 2^{K^-p^)i r(2v/^^) ^ 

T{^-q + ip)r(i + q + y^^^) 
(2Ke^)-^^ e-^, (125) 



with 



= 2 - p2 + 2p arcsin (£) , (126) 
and £ is an arbitrary non-universal constant. 

6.2 Boundary state 

Not only the vacuum amplitude, but the whole boundary state (operator) 
can be also found in the limit n — > 0. Contrary to the previous discussion 
we shall keep now the zero-mode momentum P of order 1 for n ^ 0. Hence, 
in particular, one should send p = — > in the expression H125() for 
the vacuum amplitude. With this limiting prescription X and Y sectors of 
the model are factorized completely: 

lim B = bW (g) B(y) . (127) 

Here B*^^) is the boundary state operator corresponding to the Gaussian 
boundary interaction H124|) with p = 0, and B^^^^ is the boundary state 
operator involving the y-modes only. The operator B*^^-* admits a simple 
representation in terms of the oscillatory modes Xj ()89() [41]: 

bW = g^rt^ ^ t{2k) (^^y^'^ e-^-'^-x (128) 

oo 

^""P ( ~ f arctanh(4) X^jXj ) , 
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where £y^ is some non-universal constant. The boundary state operator B*^^-* 
does not have such a simple form as H128() for an arbitrary value of q. Nev- 
ertheless its whole spectrum is known. As it was shown in Ref. [42], the 
eigenvectors of B^^^ in the Fock space of representation of the oscillatory 
modes Yj can be labeled by two increasing sequence of positive integers 

I <nl < < ■ ■ ■ < , I < < n2 < ■ ■ ■ < with > , (129) 

and the corresponding eigenvalues are given by ^: 



'•fc ~f 2^ 



where 

° "'^ r(i-,)r(i + , + .) UJ 

and i?y again is an arbitrary constant. 

A limiting behavior of the boundary state operator (|127|) as k ^ and 
K — > oo deserves to be discussed in some details. As it follows from Ea. ()13U() . 



lim B^y^ = (e^^'' + e-'^'i) go I , (131) 

where I is the identity operator. Clearly this boundary state operator cor- 
responds to the Dirichlet boundary condition for the field Y such that its 
boundary values are constrained to two points, 

YbIuv = ±^ . (132) 

This conclusion is consistent with a naive interpolation of the boundary 
constraint © to n = 0. It is also clear that the field X is subjected by the 
free boundary condition in the short distance limit: 

oo 

lim bW = S{P) exp (im ^ ) . (133) 

j=i 



^ Eg. 113011 follows from the formula (4.28) in Ref. [42], where 2pBLZ and Xblz should 
be replaced by the variables g + f- and — i respectively. 
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Now let us consider the infrared behavior of the boundary state operator 
(fT77|) . As K ^ oo, 



ip 



B«|™^r-v^ gDl e-"^^ (134) 

This equation immediately suggests that the field X obeys the Dirichlet 
boundary condition in the infrared limit: 

XB\m = ^ logr . (135) 

The infrared boundary condition for the field Y is far more interesting. As 
it follows from Ea. (|13U|) . 

If one looks at Eqs. (|134|) . (|136|) closely, several features stand out. First of 
all, it is expected a priori that the infrared behavior of p27p is controlled 
by some scale- invariant boundary state (operator). More precise, the scale 
dependence of the infrared boundary state (operator) is allowed in the form 
of a non-universal factor exp(— const k) with k = E^,R. Such factors are 
indeed presented in Eqs. (|134j) . (|136|) . The additional scale depended factor 
in (|136|) shows that when the RG "time" t = \ogn increases the field 
Y "flows" uniformly with t. It can be made into a scale invariant fixed 
point by an appropriate redefinition of the RG transformation, namely by 
supplementing it with a formal field redefinition {X,Y) {X,Y + 5t). 
This corresponds to the following modification of stress-energy tensor: 

TiR = -{dXf -{dYf + ^d^Y , (137) 

TiR = -{Bxf - {dYf + ^ a^y , 

v2 

or, in a stringy speak, to introducing a complex linear dilaton -D(X) = y. 
In a view of the non-unitary of the IPH model in the Minkowski coordinates 
(cj, ir) the appearance of complex dilaton is not particularly surprising. 

7 Exact boundary amplitude 

In this section we propose an exact expression for the boundary amplitude 
(|19j) . The expression is in terms of solutions of special second-order Ordinary 
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Differential Equation (ODE). Similar expressions are known in a number of 
integrable models of boundary interaction, beginning with the work of Dorey 
and Tateo [43] (see, e.g., [1] and references therein). Our proposal extends 
this relation to the IPH model. In this case no proof is yet available, but 
we will show in this section that the proposed expression reproduces all the 
properties of the amplitude described above. 

7.1 Differential equation 

Consider the ordinary second order differential equation 



dx^ 



p'^ + 2qKe' + K^ (e^^ + e""^) 



^'(x)=0, (138) 



where p, q are related to the components of P = (P, Q) in (|T9|) as in Eq. (|lfl6|) . 
and K = E^:R. With this identification in mind, below we always assume 
that K is real and positive. 

Let ^'-(x) be the solution of (|138j) which decays when x goes to — oo 
along the real axis, and ^'+(x) be another solution of (|138() . the one which 
decays at large positive x. We fix normalizations of these two solutions as 
follows, 

^ (2/t)~^ exp (^ - ^ e""^ ) as x ^ -oo , (139) 

and 

^ (2k)"''"5 exp I - (g + i) X - Ke"" I as x ^ +oo . (140) 

Let 

W[^'+,^'_] = ^'+ — (141) 
dx ax 

be the Wronskian of these two solutions. Then, our proposal for the function 
is 



Z = gl r-''^ , W[M/+, . (142) 

r(i -q + ip) 

To make more clear the motivations behind this proposal let us discuss some 
properties of the solutions of the differential equation H138|) 

^°Ea. l|138|l differs in tlie term 2qKe^ only from the ODE describing the boundary sinh- 
Gordon model (see Eq.(40) from the work [46] with < 0). Because of this the subsequent 
analysis is parallel to the one from the unpublished work ALB. Zamolodchikov [35]. 
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7.2 Small n expansion 



The ODE H138|) has the form of one-dimensional Schrodinger equation with 
the potential defined by the last three terms in the bracket in H138() . When 
K goes to zero this potential develops a wide plateau at 



log K ^ X ^ — log K , 



(143) 



where its value is close to —p^. In this domain each of the solutions ^-|- and 
is a combinations of two plane waves: 



*± = D±{p, q I k) e'P'' + D±{-p, q \ k) e-'^" 



(144) 



Let us consider the amplitude in (|144|) for the solution first. In 
this case it is convenient to change the variable, x = y — log(2K), and bring 
the equation (|138j) to the form 



d 



*(X) = 0; 



with 



2" k""*"^ e~^y 



Let us assume for the moment that 

- 2 < 3f?en < . 



(145) 



(146) 



(147) 



Then the term 5V^ can be considered as a small perturbation for any — oo < 
y < +00. At the zero perturbative order we dismiss SV^ in (|145jl . It gives 
an equation which can be brought to the form of Kummer's equation [24] 
by a simple change of variables. The unperturbed solution reads explicitly, 



T- e-'Py U{^ + q-ip,l-2ip,ey) . 



(148) 



The normalization of is chosen to match the asymptotic condition ()14U() . 
Now one can systematically develop the standard perturbation theory for 
the solution Therefore D+{p, q \ k) p44jl admits the following power 

series expansion 



D+ip,q\K) = {2k] 



ip 



r(-2ip) 



1 + J^d^^(p,<?) . (149) 



i=l 
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The first order perturbative coefficient d^^^ can be calculated in the closed 
form. Indeed, it is represented by the following integral: 



^1 ^^''^^ = r(i-2ip) / dzz— ^e"^x (150) 







M(i + g - ip, 1 - 2ip, z ) + q + ip,l + Tip, z) , 

which converges for — 1 < ?R.en < and can be expressed in terms of the 
generalized hypergeometric function at unity (see Eq. (7.625) in Ref. [49]): 

r(l - 2ip)r(i + q-ip-n) 
3-^2 + Q-^P, -n, -n-2\p; 1 - 2ip, ^ - n + q - ip : l). 

Evidently the power series ()149|) converges for any complex k and defines 
an entire function of for complex n ^ —2 + ^ [k = 1, 2, 3. . .) from 

the strip H147() There is no reason to expect the convergence for real 
positive n. However it is expected that for n > Ea. ()149() is valid in a sense 
of K ^ asymptotic expansion with the zero convergence radius. 

Similarly one can derive an expansion of the amplitude D- (|144|) for the 
solution In this case we make a change of variable x = - z + - log(— ); 
Eq. p38j) then takes the form, 



^{x) = 0, (152) 



with 



For 



n + 2 rv , 2{n + 2) 



6V. = ^C-^)^ e^ + f^j " ev. (153) 



4z 



3f?en<-2, (154) 

6V- can be treated as a perturbation for any real z. The small parameter 
now is . At the zero order one has 

V vrn n 

^^For q — Q the function _D+(p, 0|k) coincides up to the overall factor with a vac- 
uum eigenvalue of the Q-operator from the work [42] upon the parameter identifications: 
0BLZ = -f, Pblz = ^5 P and \%lz = -^^^j^^ In this connection, it is perti- 

nent to note that the higher-order coefficients in 11491 1 can be calculated by the method 
developed in Appendix of the work [50]. 
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where K^{x) is the conventional MacDonald function. Again, the normal- 
ization of the solution is chosen to match ()139() . The perturbative expansion 
for D_ (p, q I k) ()144|) has a form: 

D_ip,q\.) = {^y- -j0{l + Ed^\p,Q)-'-^) . (156) 



vrn 



This series converges for any complex k if n belongs to the complex half- 
plane H154|) . For real positive n it is an asymptotic series with the zero 
convergence radius. The low order coefficient d[ ^ can be easily obtained by 
direct perturbative calculations: 

d[-\p,q) = 2q(-) " ^ . (157) 



The asymptotic expansions ()149() . ()156|) are useful in study small k be- 
havior of the Wronskian (|141|) . Indeed for small k the Wronskian can be 
calculated at any point from the domain H143|) . where ^± are combinations 
of two plane waves (|144j) . Hence the Wronskian p41|) is written in terms of 
the coefficients in D± {p, q) as follows, 

W[^'+,^'_] = -2ip X (158) 
{D+{p,q\K)D^{-p,q\K) - D+{-p,q\K) D_{p,q\K)) . 



Then (|142|) leads to the short distance expansion of the partition function 
of the form ()113() . where P, Q related to p, q through Eas. (|lU6|) . The 
coefficient f-^^ in the double power series ()114() are simply expressed in 
terms of d[.^^: 

fif = d^\±p,q)d^~\TP,q), (159) 

with d^^^ = 1. In particular /-[^^(P, (5)|p=o = d^^\o,q), where d^^\p,q) 
is given by H151() . This is in a complete agreement with the result ()117l) of 
perturbative calculation in the IPH model. 



7.3 Semiclassical domain k <^ 1 <^ n 

When K is small and n ^ 1 the above expansions have to be collected. This 
regime corresponds to the semiclassical domain in the IPH model considered 
in SectionEI 
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First, let us consider the case when p and q in H138|) are order of 1; it 
corresponds to the case of "hght" insertion in H54|) with {P,Q) ~ (|61() . 
In this regime the term ^-nx equation H138() has the effect of a rigid 

wall at some point 2:0, i.e., to the right from this point, for x — xq ^ ^, 
this term is negligible, but to the left from xq it grows very fast, so that for 
X < Xq the solution is essentially zero. More precisely, when x is above 
but close to Xq the solution is approximated by a linear function, 

^^{x) ^ oq {x - Xq) , (160) 

where the position of the wall xq and the slope qq are given by 

xo = - log (^) + o(i) , (161) 



n \ n 



I Til 



as n — > 00. Here is Euler's constant. We will explain these relations 
shortly. Wronskian H141|) is easily determined. Taking, for instance, any 
point X close to the right of the wall where both equations (|148|) and (|16fl|) 
are valid, one finds it equal ao times the expression (|148|) evaluated at = 
2k e^o, or, due to equation H161|) . at e^^ = 2nA(l + o(l)), where A = (^)^. 
This leads exactly to (|67|) . 

The case of "heavy" insertion, {P,Q) ~ ^/n H72|) . can be handled simi- 
larly. The term Q-n-x g^jjj treated as a wall at xq, in the sense that 
at x—xq » ^ its effect is negligible, but it dominates at xq—x » i. To study 
the vicinity of the wall we make a change of the variable, x = ^ n 
and bring the equation H138|) to the form 1)1521) . With (|152() it is evident 
that the wall located at xo = f log(^) + ^ zq + o(^), where zq is a some 
n-independent constant. Close to the walls position, i.e., at 

|x-xo|<Cl (162) 

5V- ()153|) can be approximated by the constant 

§V_^^X + iX\ where X =(-)'' . (163) 

Therefore in the domain ()162() the solution is approximated by the Mac- 
Donald function, 

V vrn 
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with 1^0 = ^ \Jv^ — 2qn\ — (nA)^. The parameter vq here is essentiahy 
the same as in Ea.()74|l. The normahzation factor in front of K is fixed by 
matching (|164|) to the asymptotic condition ()139() . For z^o ^ and z ^ 1 
1)164(1 reduces to (|16U() . The Wronskian (|141|) can be evaluated in the domain 
^ <C a; — xo ^ 1, where both equations (|148j) and (|164jl are vahd; the result 
is exactly (|69|) with the extra factor (|78() added. 



7.4 Large k expansion 

The case of large k is expected to describe the infrared limit of the IPH 
model. For the differential equation (|138|1 it is the domain of validity of the 
WKB approximation. Applying the standard WKB iteration scheme [51] 
one finds for the Wronskian p41|) . 

logiy = lim <K / dx (V{x) — ) — ne^ — qL \ + 

where 



^(x) = ^/e^^ + e-"^ + — e^ - ^ . (166) 



nx 



The subtraction term e~~ in the integrand derives from the asymptotic 
conditions (|139j) and ensures convergence of the first integral in l|165|) as 
X — > — oo. The terms depending on L cancel the divergence of this integral 
as X — > +00. It follows from the asymptotic condition (|14()|) . Eq. ()165)) 
generates asymptotic expansion of the partition function (|19p , 

where 



Zm = gl r'^'' 2^ \. , exp ( - ^) , (168) 

and the constant kq reads explicitly, 

'^o = ^4vl ■ (169) 

^ i~2(n+2)M U - ;i+2J 
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The functions p67() are polynomials of the variables and Q of 

the degree s + 1. Their highest-order terms follow from the first integral in 

(Cnsi), 

^(norm) ^ izllV^ x (170) 

2s 

2'"(n + 2)fr(l-? + ^g^)p2/Qm 

^ n'"'r(--^) min 

which is in agreement (up to an overall normalization) with the highest-order 
terms of the vacuum eigenvalues ()95() . It is also straightforward to generate 
the full polynomials evaluating the integral ()165|) order by order in k~^. 
With the adjusted overall normalization of the local IM, this calculation 
exactly reproduces the vacuum eigenvalues 

7.5 n ^ limit 

There is no any problem with the limit n — > for the differential equation 
(|138|) and for the solution ^'-(_ (I14U() . If n = the equation turns to be of 
Kummer's type [24] and 



^+{x)\n=o = e— ^ (2Ke-)-^-^' x (171) 
U{^ + q- y^K^-p"^, 1 - 2Vk2 -p2, 2Ke^ ) . 

However, the asymptotic condition ()139() is singular as n — > and the lim- 
iting behavior of the solution is a more delicate issue. To proceed with 
the limit we will use the WKB approximation for \I'_: 

^^-^'^\x) = -^=L= exp dt{V{t) - e-f ) - f e-^ }, (172) 

where V{x) is given by (|166|) . Let us consider the argument { . . . } of the 
exponential in ()172() as n — > 0. It is easy to see that { • • • } — 
^ S{k,p), where x-independent function S{k,p) has the form (|126|) . while 
F{x) —>■ x\J — p^ as X — > +00. Thus we see that 

^'_(x) ^ ^'L'''^)(x) 2-1 (^2 _p2)-i ^ n^O, (173) 



where is a solution of the differential equation (|138|) with n = such 

that ^f?^^^^ qX^k?-p'^ as X ^ —oo. It is straightforward now to calculate 
the boundary amplitude (|142jl . The result coincides with Ea. (|125j) . 
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8 Integrable structures of the theory 



We have now seen that the expression ()142|) passed successfuhy ah available 
checks. In this section we will discuss properties of the vacuum amplitude 
1)142(1 to reveal some standard integrable structures [11,13,14] of the theory. 

8.1 Thermodynamic Bethe Ansatz equations 

It is well known from the global theory of linear ordinary differential equa- 
tions [44,45] that monodromic coefficients of linear ODE, like the Wronskian 
(|14H1 . satisfy the difference equations as functions of parameters. For the 
equation (|138() with q = and n < the corresponding difference equations 
were derived in the work [46] (see also [47,48]). The case q = 0, n > was 
studied in the unpublished paper [35]. Here we describe basic properties of 
the Wronskian (|14H) for n > and an arbitrary q. For this purpose it is 
convenient to modify slightly our notations. In particular we introduce now 
the parameter 9 such that 

where the constant kq is given by (|169|) . Also, to emphasize the dependence 
on 9 and q = ^'^"'"^ Q explicitly, we will denote the Wronskian ()141|) by 

Wq{9). Notice that there is no need to indicate the dependence on p = ^ P 
explicitly. 

The following properties of the function Wq{9) in the case n > are 
readily made using the general theory [44] and our previous analysis: 

• Wg(0) is entire function of 9. It is also entire function of the parameters 
q and p. 

• Wg(0) satisfies the so-called "quantum Wronskian" condition (see Ap- 
pendix B for details): 

W_,(0 - f )Wg(^ + f ) - Wg(^ + if^)W_,(0 - = e"'"^ ' (175) 
where a = 

• As a function of the complex variable ^, Wq{9) does not have zeroes 
in the strip j^m^] < ^ + e for some finite e > 0. 

• For \'^m9\ < f + e and 'Re9 ^ +oo: 

Wq{e) = (2^ Koe^) exp ( - e^ + 0{e-^) ) . (176) 
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For ^eO^-oo: 



ipS(n + 2) ip0(n + 2) 

Wg{e) ^ Lg{-ip) + Lgiip) e — , (177) 



where 



n _h 2h (n+2)h r(2/i)r(l + 



Lg{h) = \hr- 2"^ ™~ ('^0 



•2h\ 
n > 



Introduce the functions £±{9): 

e±i^5 {9 ± f ) {9T'f) = 1 + e""±(^) . (178) 

With the foregoing analytical conditions it is straightforward (see [52,53]) 
to transform the difference equation H175|) into a system of two integral 
equations for £±{9): 

2sin {^,) e" = ei(a) ± Si + J^^ ^ { ^^^{a - /3) x 

log (1 + e-^±(^) ) + (a - /?) log (l + e-^^C/?) ) | , (179) 

where the kernels are given by ipa-a' (o) = — i da log Sq-o-' (a) with 
sinh(f - ^) 

5'++ (a) = ^ , S-i (a) = S'++(i7r — a) . (180) 

sinh(f + ^)' 

The integral equations ()179() should be supplemented by an asymptotic con- 
dition for £±{a) as a — > — 00 which follows from Ea. p77j) . For instance, for 
Qmp < 0: 

£±{a)\a-^-oo ip ^^^^ a =F iTTg + log (Lg(ip)L_q(ip)) . (181) 

As it follows from Ea. H178() the Wronskians W±q{9) are expressed in terms 
of the solutions of (|179j) . (|18H) by means of relations: 

log W±,(0) = Tglog (2^Koe^) - + ^ ^ { 

log (1 + e-'^+(")) log (1 + e-^-(")) / 1 + e-^+(") 

^ ' ^ ^ =F i log ' 



cosh(6l-a) Vl + e-^-(°)/ cosh(6' - a) 

(182) 



39 



It is remarkable that the system ()179() differs only in the structure of "source 
terms" from the TBA system associated with the complex sinh-Gordon 
model, which is a non-compact version of the Lund-Regge model [7,54]. 
It is well to bear in mind that the classical Lund-Regge model, introduced 
in Ref. [55], is a representative of the AKNS soliton hierarchy. Notice also 
that if q = 0, then e+ = e_ and (|179() turns into the integral equation 
which describes a vacuum boundary amplitude in the boundary sinh-Gordon 
model [35]. 



8.2 T- operator 

Let us consider the differential equation in the form (|145)) and the solution 
^_ (|139|) as a function of the variable y = x + log(2K) and the parameters 
(|174|) and q, i.e., \I'_ = ^-{9,q;y). Notice that the transformations 6 
6 ± ^ leave ODE (fHS)) unchanged. Hence, the functions {0 ± ^^-.q; y) 
solve this equation as well as '^-{9,q;y), and the Wronskian 



q)=iw\^^{9+^,q;y),^.{e-^,q;y) 



(183) 



does not depends on the variable y. It can be shown (see Appendix B for 
details) that function T{6,q) is expressed through the Wronskian Wq{9) 
(HHI as 



W,{9) T{9, q)=W,{9 + ^,)+W,{9-^) , 



(184) 



and 



T{9 + i7r^,q)=T{9,-q) . (185) 

Since the vacuum amplitude ()142|) differs by ^-independent factor from the 
Wronskian^^, one can replace Wq by Z in (|184|) . It is also clear that the 
function T{9,q) admits 9 — > +oo asymptotic expansion similar to (|167|) . 
The foregoing in turn leads us to 



CONJECTURE 



There exist an operator T which acts invariantly in the Fock space J-p and 
satisfies the following conditions: 

^^Here we assume that the parameter r in l)142|l does not depends on 9. This assumption 
follows from the normalization condition I65II . 
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9(n + 2) 

Being considered as a function of the parameter A ~ e " 

T = T(A) admits the convergent power series expansion of the form: 

oo 

T(A) = 2cosh(^)+^GfcA^ (186) 

fc=i 

The coefficients in H186|) satisfy the condition, 

MGfcM = (-1)'' Gfc , (187) 

where the operator M : ^(p^q) J^(p._q) flips the overah sign of the 
field dY, i.e., MaFM = -dY. 

The operator T(A) commutes with all the IM from AKNS series and, 
hence, 

[T(A'), B(A)] = . (188) 

Here the boundary state operator B is understood as a multi-valued 
function of the spectral parameter (|16j) . The eigenvalue of T corre- 
sponding to the Fock vacuum | P ) G J-^p coincides with T{9, q) H183|) . 

The boundary state operator B and T satisfy the T—Q Baxter equation 

The operator T(A) is a generating function of the AKNS series of local 
IM. In more exact terms the local IM are generated in the large A 
asymptotic series expansion of T(A): 

T(A) = exp - 27rz^(Ae'^) ) +exp (^27rz^(Ae"v) ^ , (189) 

where 

V^r(^) ^ ' 

The local IM li"""™^ in (ITO are normalized in accordance with the 
condition. 



1 / 1 I ^ V 2 ' n+2 ) 



n 2 ^ Jo I 2l+m=s+l nU) 

\f-i {dX)^{dYy 
n J l\ ml 



X 
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where omitted terms contain higher derivatives of 5X 



Eventually, the operator T(A) should be viewed as a quantum version of the 
transfer-matrix T(A) ()5U|) for the AKNS linear problem (|46l) with j = \- 
Following the line of Ref. [11] it seems possible to express the coefficients 
in the power series expansion H186() in terms of 2A:-fold integrals over chiral 
vertex operators involving the holomorphic component of field X. Since the 
nonlocal operators Gk commute among themselves [Gfc,Gm] = and also 
commute with all the local IM I^, [Gfc,!^] = 0, they are called the nonlocal 
Integrals of Motion. Unfortunately the "explicit" formulae for nonlocal IM 
are not particular useful neither to prove the above listed properties of the 
operator T nor for calculations of its spectrum. For this reason we do not 
present these formulae here. 

Finally we note that the vacuum eigenvalues C^™'^^ of nonlocal IM, 

GfclP) =g1""'^|P) , (192) 

can be algebraically expressed in terms of the coefficients ^ (p, q) in the 
formal power series expansion For instance, 

G^r^ = -4 sin (^) sin d^^\p,q)= (193) 

_!i£(i+i) I 

^/^r(i + i) r(i + i-f^)r(i + i + ^) ' 

8.3 Commuting families in the quantum AKNS hierarchy 

In Refs. [11,56] the quantization procedure for the KdV hierarchy was devel- 
oped. In the BLZ approach quantum transfer-matrixes are defined in terms 
of certain monodromy matrixes associated with 2j + 1 dimensional represen- 
tations of the quantum algebra [/q(sZ(2)). The similar C/q(,s?(2))-structure 
can be observed in the quantum AKNS hierarchy. In consequence of this the 
quantum KdV and AKNS hierarchies share some common formal algebraic 
properties. In particular, the quantum transfer-matrixes Tj associated with 
the 2j + 1 dimensional representations of f/q(s/(2)) in both hierarchies are 
recursively expressed through the operators T = Ti and the unit operator 

2 

I = To by means of the same fusion relation (see, e.g., [56]): 

T(A)T,-(q^-+h)= T 1 (q^+^A) + T^.^i (q^) . (194) 

""^^Note that the local IM Ij"""^™-* are defined unambiguously with the normalization 
condition H191|l . 
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Here q is the parameter of deformation of C/q(s/(2)). For the AKNS hierar- 
chy it should be chosen as 

q = e"^ . (195) 

In Section lST^ the vacuum eigenvalue of T was identified with the Wronskian 
H183|) . With the fusion relation H194|) we can express now the all vacuum 
eigenvalues, 

Ty(A)|P) =rj™')(A)|P) , (196) 

in terms of solutions of the differential equation H145|) . The result appears 
to be in a remarkable form generalizing (|183j) : 

(A) = i(-l)2^+iW [e + ^^iggi) , q- y) ,^.{e- ^iggil ,q-y)] .(197) 

Recall that A = (^ e^)~^, where the constant kq is given by (|169|) . 

We can not resist the temptation to mention here an evidence of existence 
of additional commuting family in the quantum AKNS hierarchy which does 
not have a classical counterpart. Indeed, let us look at the solution ^'-|_ (|14U() . 
On the complete analogy with the discussion from previous subsection we 
consider ^'_|_ as a function of the variable y = x + log(2K) and the parameters 
9 (|174j) and i.e., \I'+ = ^+{9, q; y). It is easy to see that the transformation 

y^y±i7,, e^9±^, q^-q (198) 

leaves the equation H145|) unchanged. Because of this the functions ^'-i- (^ it 
^q^, —q; y ± ivr) solve (|145j) and the following Wronskian does not depend 
on y: 

f{9, g) = i W [ ^+ (0 - -g; y - ivr) , (0 + ^, -g; y + ivr) ] . (199) 

It is straightforward to show that T(9,q) satisfies the relations: 

Wgi9) fi9, q) = e-^^5 W.,{9 + ^) + e^"" W.,{9 - j^) , (200) 

and 

f{9-^,q)=f{9,q) . (201) 

The evident similarity between (|MH) . !^UT^ and (fTH^ . ((TH1|) suggests to in- 
terpret T{9, q) as a vacuum eigenvalue of "dual" transfer-matrix T. Due to 
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Ea. ()2Ul() . it is expected that T admits a convergent power series expansion 
in terms of the "dual" spectral parameter A = k""^^. Using properties of 
the Wronskian ()199|) . it is easy to guess a set of general conditions for the 
operator T(A) which is similar to the one expounded in Section lH^ for T(A). 
In all likelihood the appearance of "dual" transfer-matrix is related to the 
hidden [/q(si(2|l))-structure in the quantum AKNS hirarchy^^. 

9 Infrared fixed point of the IPH boundary fiow 

In this work we have mainly focused our attention on the vacuum boundary 
amplitude ((T^. Of course, even the exact expression H142|) does not de- 
fine the boundary state completely. However, as we saw in Section|Hl some 
properties of the vacuum eigenvalue inherit important general features of 
the whole boundary state operator. Among them the remarkable T — Q 
equation which is a keystone relation in the theory of integrable systems. 
Another property of Z which should be understood in more general terms 
is the large k expansion (jlOTj) . It suggests that the boundary state oper- 
ator B not only commutes with the local IM, but it admits an asymptotic 
expansion in terms of these operators, 

I' oo -jj{norm) 

B.B,Kexp|-|:5-l^(^)'|, (202) 

where the local IM 1^°°'^™'* are normalized in accordance with the condition 

(Unii)- 

From the physical point of view the most interesting object to be dis- 
cussed is the infrared boundary state associated with the operator Bjr pfl2j) . 
According to our consideration this boundary state should have the form 

|5)iR = ^d2p |Tp) , (203) 

where the amplitude Zjr is given by Eq. (|168() and therefore the states | Xp ) G 
J'P <^ are normalized as follows: 

( P' I Tp ) = S(P' - P) . (204) 

Of course, the boundary state (|2fl3j) should obey the integrability condition 
(jEJ. It is also expected to be a conformal boundary state which is some 

^*The screening operators 18111 . 1991 can be interpreted as generators of the Borel sub- 
algebra of the quantum afhne superalgebra f/q(s/(2|l)) with q = e""^". 
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deformation to the domain n > of the infrared boundary state with n = 
described in Section 16.21 In Appendix C it is shown that a boundary state 
satisfying the above mentioned conditions must also possess an extended 
conformal symmetry generated by the holomorphic spin-1 and spin-2 cur- 
rents: 



[ W^'^^ (r) - W^'^^ (r) ] I i? )iR = = 1, 2) 



(205) 



where 



-{dY*f + — s^y* 

v2 



In Ea. ()2U6() we use the notations 



y* 



(IT) \ 

Notice that the spin-2 current W2 produces the Virasoro algebra with 
center charge c = — 2. If 



2 V 2 



(206) 



(207) 



(208) 



the building blocks (Ishibashi states) | Tp ) in (|2U3() are defined by Eas. (|2U4() . 
H205|) uniquely for the Fock spaces J-^ J--p with the zero-mode momentum 
P = (P,Q) [3,4]. In particular, 



|Tp) = exp(X;f^-fc^-fe 



k=l 
3 



1 + -^ — \ 2y*iy*i + 



( {Y*-if + ^ Yi, ) ( {Y*_,f + ^ y?2 ) + 



+ 



★ \2 



j=,Y^2)mi?-^,Y*^2) + 



V2 



(209) 



Here X'^^, Y*^ are the oscillatory modes of the fields H207|) . 

It can be shown that the integrand in Ea. H203|) . where | 2p ) given by 
H209|) , is well defined even for the exceptional values q* = it ^ , ± | • • • and 
therefore, we propose ()2U3() . H209|) for the infrared boundary state of the IPH 
model. Notice that it splits into the Dirichlet boundary state corresponding 
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to the Dirichlet boundary condition = const, and the n-independent 
conformal boundary state for the field 1"*. Unfortunately a consistent La- 
grangian description of the y*-component of | -B)ir is still a question for 
the authors. 

We would like to conclude the paper by a remark about the higher level 
boundary amplitudes Ba(P) in Ea. H13() . It seems likely that they also admit 
a description in terms of ODEs similar to (|lH8j) . Such higher level differential 
equations can be constructed along the line of Ref. [57] exploring analytical 
properties of the boundary state operator B as a function of variable 6 p74|) . 
We intend to return to this problem in a future work. 
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10 Appendix A 

Here we present an explicit form of the first local IM of the AKNS series as 
an operators acting in the Fock space. Let us introduce the notations, 

K:::t= E (sir ■■■{smT- -.x^i ...x^- ■. . (210) 

si+...sm=0 

The normal ordering : : in this formula means that the operators X] ()89|) 
with the bigger s are placed to the right. The local IM can be written as 
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follows: 

11 = [Xll + Xll + I[^^^\P,Q)], (211) 

12 = R-^ [^Xlll + 2nXlll + {3n + 2)QXll + nQXll + 

2n P Xll - 2i (n + 1) V^X\l + it"^ {P, Q) ] , 

and 

I3 = [nXjJii + (5n + 4)X2222 + 6nXiJ22 + 2PnXiJJ + 

2(5n + 4) Q Xlll + 6nP XHI + 6nQ Xjjg - 6i ^/^ (n + 1) X\ll + 
§ (3P2 + 3Q^ - 1) XjJ + i (3(5n + 4)^^ + SnP^ - 3n - 2) Xgg + 
6raPQ Xll - 6i (n + 1) Q X^^ - 

(n^ + 3n + 1) X}} - (n^ + 4n + 2) X?^ + 4™^) (P, Q) ] . (212) 
Here/i™^)(P,Q) are the vacuum eigenvalues of the AKNS integrals (pi)) . 



11 Appendix B 

Here we prove the quantum Wronskian condition (|175|) and Eq. 1)184(1 . 

We start with an observation that the following transformations of the 
variables {y, n, q), 

A : y—>-y + 'nT, k —>■ e"+'^ k, q ^ —q , 

n : y^y, k ^ , q^q, (213) 

leave the ODE l|145p unchanged while acting nontrivially on its solutions. 
The transformation A applied to the solution (|14U|1 yields another solu- 
tion, and the pair of functions 

= ^+{y, e, q) , = ^+{y + ivr, 9 + -q) , (214) 

with 9 given by (|174() . forms a basis in the space of solutions of ()145p . It is 
not difficult to check that, 

Ty[^+,A^+] = e'^(''-5) , (215) 

i.e., the solutions (|214|) are indeed linearly independent. The solution 
(|139|1 can always be expanded in this basis, in particular 

^_=a + 6 A1'+ . (216) 
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Using Ea. (|215p we conclude that 

b = e'^(i-«) Wg{e) , (217) 

where Wq{0) = W['^>-^-,'^^~\. The transformation A leaves the solution 
unchanged^^, i.e., ^-{y,9,q) = ^_ (y + ivr, + -g) . This allows one to 
express the coefficient a in 1)216^ in terms of the Wronskian Wq{9) as well: 

a = -ei-(|-<?) W^g{e + ^) . (218) 

Let us apply now the transformation 0, ()213|) to the both sides of Ea. (|216() . 
It is apparent that ^'_|_ is invariant with respect to the action of Cl, and 
hence, 

^ gi-(i-g) (^y4/^(0 + |izL) A^+-VF_g(6' + i7r) ^'+) . (219) 

The quantum Wronskian condition H175() follows immediately from ()216() . 
(EHI), (PT3|) and the relation 

W[if-,n'^-] =i . (220) 

It is easy to derive now from the last equation that three solutions 
r2\I'_ and r2~^^'_ satisfy the relation, 

T{e, q) = + , (221) 

with function T{9,q) defined by Eq. (|183|) . Taking the Wronskian from the 
both sides of this equation with the solution ^'^ = = il^^'I'-i- we arrive 
to Eq. (tTMl) . 

12 Appendix C 

Here we show that the infrared boundary state of the IPH model, | -B)ir, 
satisfies the conditions (|^n ^ - PT7|l . 

Our analyses is based on the assumption that | -B)ir possesses the con- 
formal symmetry, i.e.: 

[T(i«')(r)-t (i^) (r) ] I i? )iR = , (222) 

^^This property of "if- leads immediately to Eg. 118511 . 
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where the holomorphic and antiholomorphic components of infrared stress- 
energy tensor have the most general admissible form: 

2^(iR) ^ • 5X + i p • S^X , (223) 

^(iR) ^ -ax-ax + ip-a^x . 

The constant vector p = {px,Py) is unknown a priori, but it is expected to 
be a function of coupling constant n. A comparison between ()223|) and the 
infrared stress-energy tensor for n = (|137|) shows that 

Py\n=o = . (224) 

The local fields (f^^ and the first nontrivial local IM I2 (S2) (PTTj) allow 
one to construct holomorphic spin-3 and antiholomorphic spin-(-3) currents 
through the relations: 

dW^^^^ = [T,h] , BW^^^^ = [ r , I2 ] . (225) 

Since the infrared boundary state satisfies both Ea. H222|) and the integrabil- 
ity condition it is apparent that 

[ (r) - M'""^ (r) ] I S )iR = . (226) 

An explicit form of W3 is not particularly important for the present dis- 
cussion. It is important that this field and T^^^^ generate a spin-1 current 
through their operator product expansion: 

/Tr> \ 

Explicit calculations shows that has the form 



Vrf = adX + i(3dY , (228) 



with 



a = 2 py^/n {I + n + ip^y/n) , (229) 
(5 = np^ + p^(2 + 3n) + 2ip^ (1 + n)Vn -2n- 1 . 

As it follows from Eas. ()229() the current W^^^^^ may vanish for py = ot py = 
±V^ + 2 only. Since none of these is consistent with Eq. (|224j) . the infrared 
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boundary state should satisfy the condition ()2U5p for the nonvanishing spin-1 
current given by Eas- I^^ and (f^^ . 

Let us consider now the operator product expansion Wi^^\u)W^^^\v) . 
It contains a singular term ~ (n — v)~'^ which involves the spin-2 current of 
the form, 

W^^^^ = 2inp{dXf + AnadXdY + i/3(6n + 4){dYf + -^d^X + 6d^Y,{230) 
where 

7 = n{Pp:c-iapy) + i{n + l)P^/n , (231) 
5 = -a^/n{n + l + ipx^/n) + P py{3n + 2) . 



With Eas. (|228|) . H23U|) it is easy to see that the operator product expansion 
W^^^\u)W2^^\v) produces one more spin-1 current 

Wi^^'^ = 2m ap dX + { na^ - (3n + 2)P^) dY . (232) 



This current does not vanish provided exists. Therefore one needs 

to explore two possibilities: either and are linear independent 

or linear dependent currents. The first possibility implies that | B )ir is 
the ra-independent Dirichlet boundary state associated with the Dirichlet 
boundary condition = const. It should be apparently ignored. Hence 
Wj;^^^ ~ lyf With Ens. (l228|l . ([232|l it gives the relation 



Q ± Vn + 2 P = . (233) 

Now it is convenient to introduce a new set {X*, Y*) related to the basic 
fields {X, Y) through the (complex) orthogonal transformation: 



Here the sign factors are dictated by the choice of sign in 1)233^ . Since 

/TO -J 

W^'"^' ~ dX* one can choose (without loss of generality) the infrared stress- 
energy tensor (|223)1 in the form 

T(^^) = -{dX^f + W^^^^ , (235) 

where 

W^l^^ = -{dY^f + ^ d'^Y* . (236) 
v2 
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It follows immediately from Eas- lf^^ and (^H^ that cr^ = 1. 

Thus there are four nonvanishing spin-2 currents, r(ii^), T^f^), 
and d'^X*, satisfying the condition {W2 — W2) \ = 0. If one admits 

now their linear independence, then d'^Y* can be expressed in terms of 
these currents and, consequently, should also satisfy the above condition. 
This will eventually lead us to the Dirichlet boundary state corresponding 
to = const. Therefore we reject this possibility and accept that the 
current H^g^^'' is linearly expressed in terms of T^^^\ VFg^^'' and d'^X*. It is 
straightforward to check that this is indeed the case provided 

a = +l . (237) 

The remaining ambiguity in sign factors in Eas. ()234() should be resolved in 
accordance with the condition 1)2241) . In this way one arrives to Eas. H2U5() - 
(EOZl). 

Finally we note that the local IM I2 ^211^ can be expressed in terms of 
the fields (IM . (PUTI) as follows 

I2 = (2n + l)y^4"=°)[y*]-2^/2^JI2[VFf^^ VFi^""^] , (238) 
where 

t='\Y^] = ^ fj^'^idy-f. (239) 

and 

HWi, W2]= — ( (n + 1) W2W1 + ^ {Wif ) . (240) 

Jo 

Evidently (J2 — JI2) I = 0, therefore the integrability condition (I2 — 

S-2) I ^ )lR = for the boundary state (PUS)) . (PU^ with n > follows from 
the fact that it holds for n = 0. 
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